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ABSTRACT 

We investigate the evolution of the galaxy merger rate predicted by two semi-analytical 
galaxy formation models implemented on the Millennium Simulation of dark matter structure 
growth. The fraction of merging galaxy pairs at each time-step of the simulation is derived 
from the galaxy catalogues obtained by the models and the results are compared with various 
observational estimates of merger fractions taken from the literature. We find a good match 
between the pair fractions derived from the simulation and the observed counting of galaxy 
pairs obtained by different sources in the redshift range < z < 1.2. The predicted evolution 
of the number of galaxy mergers per Gyr grows with redshift as an exponential rate given 
by T mrg oc (1 + z) m , with m ranging from 0.6 to 0.8 for < z < 2, depending on the 
luminosity and mass ratios of the merging galaxies. Our results are in agreement with recent 
observational results that argue for a fiat evolution in the fraction of galaxy mergers since 
z ~ 1.2. The weak evolution predicted for the galaxy merger rate in an hierarchical model 
universe shows that the mass assembly evolution of galaxies through mergers does not follow 
the rapid evolution of the halo merger rate obtained in previous studies. 

Key words: galaxies: evolution — galaxies: formation — galaxies: interactions — method: 
data analysis 



1 INTRODUCTION 

The role played by galaxy mergers is of paramount importance 
for galax y formation and evolution . Since the seminal numerical 
work by iToomre & Toomrej dl972f) . galaxy mergers are thought 
to be relevant mechanisms in the build up of the spheroidal 
component of galaxies. Recent numerical works also argue that 
gas-ri c h mergers are able to form disk galaxies dRobertson et al.l 
|2004 iGovernato etafl 120071) . In the local Universe, observa- 
tions of remnants of galaxy mergers strengthen their effects on 
galaxy morphologies and on triggering star formation in merg- 
ing galaxies, unveiling the i r power on driving galaxy evolution 
(e.g. ljog & Mavbhatell2006l: ISchweizer & Seitzerl|2007|) . Actually, 
in a recent work lMateus. Jimenez. & Gaztafiagal J2008l) show that 
galaxy mergers are very correlated to the star formation history of 
galaxies, driving the environmental dependence in galaxy forma- 
tion by regulating the star formation process. 

The importance of galaxy mergers is highlighted when we 
consider the hierarchical build up of cold dark matter (CDM) struc- 
tures according to the standard ACDM cosmology, which forms the 
basis of the subsequent processes related to galaxy formation and 
evolution. In this picture, the structure growth is driven by the hier- 
archical accretion and merging of dark matter haloes, with galaxies 
forming inside these structures following the radiative cooling of 
bary ons and the condens ation of the gas through dissipative cool- 
ing l lWhite & Reeslll978h . It follows that the rate at which galaxy 
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mergers evolve is a direct consequence of the hierarchical growth of 
structures in any galaxy formation model based on this framework. 
Predicting this evolution from the models and comparing it to ob- 
servations of galaxy mergers at different lookback times represent, 
therefore, a key test to determine whether and how galaxies form 
by mergers or whether the models give us the correct predictions. 

On the observational side, the galaxy merger counting is a 
challenging effort. Two main approaches are used to measure the 
evolution of galaxy mergers at low and high-redshift. The first one 
is based on the counting of close galaxy pairs, which provides re- 
liable estimates of major m e rgers out to z ~ 1 



19971: iLe Fevre et al 



20001 : iBundv et alj |2004 iLin et all 2004; 



Kartaltepe et alj |2007), despite of the selection ef fects and biases 
inherent to this method (e.g. IPatton et"af I Eoooh . Other way to 
counting galaxy mergers is through detecting direct signatures of 
galaxy interactions, as visual disturbed mo rphologies or asymme- 
tries detected by automated procedures (e.g. lConselice et al]|2003l : 
|Pe Propris et alj2007l : lLotz et alj200j) . 

In this Letter, we examine the galaxy merger rate evolution 
predicted by galaxy formation models and its comparison with 
observational measurements of merger rates taken from the liter- 
ature. We use the public dat a available by the Millennium Sim- 
ulation JSpringeletaT1l2005l) to derive galaxy merger properties 
predicted by two semi-analytical galaxy f ormation models im ple- 
mented on the simulation a nd described bv lBower et al. and 
|Pe Lucia & Blaizotl ( [20071) . 

We present a brief description of the Millennium Simulation 
and the semi-analytical models used in this work in Section [2] In 
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Section [3] we describe the methodology used to estimate galaxy 
mergers from the simulation outputs. The main results of this work 
are presented in Section|4]and our conclusions are summarized on 
Section|5] 



2 THE MODEL UNIVERSES 

In this section we briefly describe the model universes used in this 
work based on one of the largest cosmological simulation of struc- 
ture growth carried out so far, the Millennium Simulation, and on 
two state-of-the-art semi-analytical models of galaxy formation de- 
veloped by independent groups. 

The Millennium Simulation of dark matter structure growth, 
recently carried out by the Virgo Consortium, follows the evolu- 
tion of N = 2160 3 particles of mass 8.6 x 10 s h' 1 M , within 
a co-moving box of size SOOfe" 1 Mpc on a side, from redshift 
z — 127 to the present l lSpringel et ai1l2005l) . The adopted cos- 
mological model is the concordance ACDM model with parame- 
ters Q m = 0.25, Q b = 0.045, h = 0.73, f2 A = 0.75, n = 1 
and er 8 = 0.9, where the Hubble constant is parametrized as 
Ho — lOO/i km s _1 Mpc -1 . Haloes and subhaloes in 64 output 
sn apshots were identifie d using the SUB FIND algorithm described 
in ISpringel et all fcOOll) . and merger trees were then constructed 
that describe how haloes grow as the model universe evolves. These 
merger trees form the basic input needed by the semi-analytic mod- 
els described in the following. 

We explore the galaxy catalogues obtained by two semi- 
analytical galaxy formation models implemented on the Millen- 
nium Simulation. We use the results obtained by th e model pre- 
sented and discussed by |Pe Lucia & Blaizol ||2007|), which is a 
modified version of the models used in Croton et al. I j2006h and 
|Pe Lucia et alj J2006l) . The galaxy catalogues generated by this 
model have been made public on the Millennium download site 
at the German Virtu al Observatorfl We also use the model de- 
scribed in detail by iBower et alj d2006l) . It is an extensio n of 
the GALFORM m odel implemented by ICole et alj ( feOOOl) and 
iBenson et alj d2003l) . The catalogues generated by this model are 
publicly available on the Millennium download site at the Univer- 
sity of Durham^ 

These semi-analytical galaxy formation models have signif- 
icant differences in the general procedure to form galaxies from 
the dark matter haloes. First of all, many aspects of the baryonic 
physics involved in the galaxy formation process are treated in a 
completely different manner, mainly those related to the energy 
feedback from active galaxy nuclei to regulate the star formation 
in galaxies. In addition, the construction of the merger trees from 
the dark matter halo distribution, which form the basis of the semi- 
analytical approach integrated on the numerical simulations, has 
been done in independent ways. 



Bower et al. d2006h model and about 4.35 x 10 6 galaxies from 



De Lucia & Blaizotl l [2007h model. For all galaxies in these cata- 



3 GALAXY MERGER RATES 

We derive the predicted evolution of the galaxy merger rate from 
the semi-analytical galaxy formation models implemented on the 
Millennium Simulation. We start selecting all z — galaxies 
with rest-frame absolute magnitudes Mb ^ —18, resulting in 
two base catalogues containing about 3.33 x 10 6 galaxies from 
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logues, we recovered their complete merger trees by using the post- 
processing facilities carried out on the Millennium databases. 

The first step to derive the model predictions for the frac- 
tion of galaxies undergoing mergers is to identify the progenitors 
of each galaxy selected from the merger trees at each time-step t 
of the simulation. We then define that a galaxy has experienced a 
merger if it has at least one pair of progenitors at the previous time- 
step. In other words, we are able to select galaxy pairs at a given 
t that will merge to form a single galaxy in the immediately pos- 
terior time-step. This description is particularly interesting if one 
wants to compare the results from the simulation with observational 
estimates of merger fractions bas e d on pair count statist ics (e.g. 
IPatton et alj|2002l ; iLin et"ai1l2004k [Kartaltepe et af]|2007l) . In our 
analysis, the progenitor pair of a merge-product galaxy is composed 
by its main progenitor (defined as the progenitor with the largest 
stellar mass), selected to be brighter than a B-band absolute mag- 
nitude limit, Mg m , and an ordinary progenitor with Mb ^ —18. 
In our analysis, we considered two magnitude limits for the main 
progenitors, M]J m = -19.5 and M l ^ m = -20, in order to in- 
spect how this choice affects the derived pair fractions. We have 
also adopted an additional approach to select galaxy pairs by con- 
sidering only the ordinary progenitors with stellar masses larger 
than 1/4, 1/6, and 1/10 of the main progenitor stellar mass, in- 
dependently of their luminosity. In such analysis, we selected only 
the main progenitors with Mb < —20. This procedure allows us 
to investigate how the estimated pair fractions depend on the mass 
ratio of the pair components. 

We define the pair fraction as the number of galaxies in the 
progenitor pairs (main progenitors plus companions), N p (z), di- 
vided by the total number of main progenitors brighter than M^'" 1 , 
Nt(z), at each time-step of the simulation: f p (z) — N p (z)/Nt(z). 
The merge fraction is directly derived from the pair fraction by con- 
sidering only the total number of progenitors pairs that will merge, 
that is, fmrg(z) — 0.5/ p (z). We can directly relate this predicted 
merger fraction to observations of close galaxy pairs which are fre- 
quently used to derive observed merger fractions. From observa- 
tions of physically close galaxy pairs one can determine the number 
of close companions per galaxy, N c , which is identical to the pair 
fraction in a volume-limited sample with no t riples or h i gher o rder 
N-tuples. Following the expression adopted bv lLin et"aT1 ( l2004l) . the 
galaxy merger rate per Gyr estimated through pair count statistics 
is given by 



Rrarg — (0.5 + G) N c (^z) P T nrgT rnrg 



(1) 



where N c (z) is the average number of companions for galaxies 
within the observed magnitude range, P mrs denotes the probabil- 
ity that galaxies in close pairs will merge (typically assumed to 
be 0.5), and T mrg is the time-scale for physically associated pairs 
to merge. The factor 0.5 is to convert the number of galaxies into 
the number of merger events (assuming only pairs, ignoring triples 
and higher order N-tuples). The physical galaxy pairs are expected 
to merge on the time-scale T mrg which depends on the intrinsic 
characteristics of the pairs, such as the relative mass ratio, orbit pa- 
rameters, and on th e structure of the merging galaxies, thus varying 
for each pair (e.g. ICon selice 2006 1. In general, one assumes that 
T mr g ~ 0.5 Gyr, a value suggested by simplified models based 
on the dynamica l friction time-scale and also predicted by N-bod y 
simulations (e.g. lMihos|[l995l : |Patton et al.ll200d:IConselicdl2lj06h . 
In equation {Q we also inc lude an additional parameter, G, in- 
troduced by |Linetal.N2004l) to account for differences in sample 
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selections when comparing merger rates estimated by pair counts 
and by morphological ap proaches. Here w e will adopt the value of 
G = 1.24 determined bv llin et all d2004l) . 

In the case of the simulation, a merger between one or more 
ordinary progenitors and the main progenitor results in the for- 
mation of a single galaxy at a posterior time-step. Thus, we can 
evaluate a galaxy merger rate (number of galaxy mergers per Gyr), 
Vm r g(z), by dividing the merger fraction at each z by the time in- 
terval between the two time-steps associated to the progenitors and 
to the new formed galaxy: r mrg (z) = fmr g {z)St~ 1 . In the red- 
shift range < z < 1, the mean value of St is 0.34 Gyr, implying 
that the progenitor pairs of a given galaxy can be considered as 
physical pairs. Thus, the predicted galaxy merger rate can be easily 
compared with that derived from observational data, R mrg , based 
on pair count statistics. 



4 RESULTS 

In Fig.[T] we show the predicted evolution of pair fractions derived 
from galaxy formation models, f p (z), and a compilation of obser- 
vational results in the redshift range < z < 1.2. For the simu- 
lation, we show results of pair fractions obtained assuming mag- 
nitude limits to select the main progenitors of Alg m = —19.5 



0.10 



(bottom lines) and M™ 



-20 (top lines). The observational 



data points shown in this fi gure were obtained f rom estimates us - 
ing galaxy pair count s by IPatton et alj d2002h . iLin et al. I J2004h . 
iBell et all J2003) . and |Pe Propris et al.1 J2007I) . For z < 0.3, we 
show the close pair fractions obtained by Patton et al. from counting 
close dynamical pairs of the CNOC2 survey, with luminosities in 
the range —22 < Mb < —19, and the result obtained by De Pro- 
pris et al. for a galaxy sample with —21 < Mb — 5 log h < — 18 
and 0.010 < z < 0.123 drawn from the Millennium Galaxy Cat- 
alogue. At higher redshifts, we show the mean value of the close 
pair fraction obtained by Bell et al. after analysing 0.4 < z < 0.8 
galaxies with M B < -20 from the COMBO-17 survey and the 
pair fractions obtained by Lin et al. from analysis of galaxies with 
luminosities -22 < Mb < -20 and redshift range 0.5 < z < 1.2 
from the DEEP2 redshift survey. It is worth noting that we have ad- 
justed the values of Lin et al. downward by 65 per cent to account 
for pr ojection within galaxy groups, as suggested by IBell et all 
d2006h . 

There is a good match between the e volution of the pair frac- 
tions predicted by the lBower et alj | |2006|) model and that inferred 
by the observational data, considering the mag nitude limits adopted 
to sele ct the main progenitors. The model bv lDe Lucia & BlaizoJ 
d2007l) underpredicts the pair fractions at z > 0.5, showing a de- 
creasing evolution which is clearly inconsistent with the behaviour 
expected from the observations. The evolution of the pair fractions 
derived from the Bower et al. model is strong at z < 0.3, and be- 
comes weak at higher redshifts, independent of the value of M^" 1 
adopted to select galaxy pairs. The amplitude of the pair fractions is 
highly dependent on the value of Mg m , being higher for pairs con- 
taining brighter galaxies. This trend is associated with the strong 
dependence of the galaxy growth throug h mergers on stellar mass 
dMaller et alJl2006l : iGuo & Whitell2008l) . We have also compared 
our resul ts with the recent measu rements of close pair fractions ob- 
tained bv lKartaltepe et all d2007l) . who claims for a strong evolution 
of the pair fractions since z ~ 1.2 after analysing the galaxy popu- 
lation in the COSMOS field. They adopted photometric redshifts to 
select galaxy pairs in the redshift range 0.1 < z < 1.2, with a mag- 
nitude limit Mv < —19.8. We found a good agreement in the pair 
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Figure 1. Evolution of the pair fractions predicted by the semi- 
analytical galaxy fo r matio n models of iBower et alj j2006l solid lines) and 
|Pe Luci a & Blaizot 1 2007, dashed lines) and the observational data taken 
from different sources for galaxies in the redshift range < z < 1.2. The 
lines correspond to model predictions considering Mg m = —19.5 (bot- 
tom lines) and Mg m = — 20 (top lines) to select the main progenitors, as 
described in Section[5] 



fractions for z < 1. For higher redshifts, there is an excess in the 
pair fractions derived bv lKartaltepe et all d2007l) . inconsistent with 
the weak evolution of f p (z) predicted by the Bower et al. model. 
This discrepancy can be due to the large errors in the photometric 
redshifts used to select the pairs and the consequent projection ef- 
fects, which become more important at higher redshifts. We hope 
that a future spectroscopic follow-up of their data and other pair 
fraction estimates from high-jz surveys will help to clarify this is- 
sue. 

In Fig. [2] we show the evolution of the merger rates predicted 
by the models, r mrg (z), and the merger rates derived from the 
observations. Galaxy merger rates are estimated from the obser- 
vational data considering the time-scale T mrg = 0.5 Gyr for the 
measurements from pair counts, assuming by simplicity that this 
value is characteri stic of the merging time-scale for close pairs (e.g. 
IPatton etal]|2000l) . as discussed in Section [3] For this compari- 
son, we also added the merger rates estimated from merger anal- 
ysis based on galaxies showing morphological distortions. This ap- 
proach is significantly distinct to the pair counts a nalysis, since the 
merg ers are observed after their occurrence (see |Pe Propris et"al] 
1 20071 for a discussi on between thes e two methods). We show the 
results obtained by IConselice et all ( 120031) from a morphological 
analysis of galaxies in the Hubble Deep Field. They have detected 
ongoing major mergers by measuring an asymmetry parameter 
from galaxy images. We also show the recent measu rements of 
morph ologically identified galaxy mergers obtained by iLotz et all 
d2008l) from the analysis of galaxies brighter than Mb = —20.5 
in the HST survey of the Extended Groth Strip. Here we adopted a 
time-scale of 1 Gyr to derive the merger rates from these morpho- 
logical studies. 

As in Fig. [T] there is a good agreement between the observa- 
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Figure 2. Evolution of the galaxy merger rate p r edicte d by the semi- 
analytical galaxy fo r matio n models of iBower et al solid line) and 

|Pe Lucia & Blaizol dashed line) and the observational data taken 

from different sources for galaxies in the redshift range < z < 2 (see 
Section|4]for detailed description). The lines are as in Fig.[T] 



tional data points and the Bower et al. model predictions, consider- 
ing the error bars of the measurements and all the uncertainties in- 
volved in deriving the observed merger rates (e.g. time-scales). The 
evolution of merger fractions derived from observational data sets 
is usually fitted with a function / oc (l+z) m , where the value of m 
varies from about to 4 (e.g. iLe Fevre et al . 2000; C arlberg et alj 



120001 : IConselice et"ai]|2003l: iLin et alj|2004 Lotz et alj|2008h . In 
our analysis, the fitting of the Bower et al. model predictions with 
the same function, considering only the points in the redshift range 
< z < 2, gives: 



T mrg {z) = 0.029(1 + z) - 767 for M^ m = -19.5, and 



mrg \ 

r mrg (z) = 0.046(1 + z) ' 596 



for Mg m = -20.0. 



(2) 



The m ~ 0.6 — 0.8 value obtained for this model is consis- 
tent with recent results obtained by different approaches, which 
gives support for roughly a flat evolution of the merger rate 
bundv et alj|2004 ILin et al ll2004 Eotz et alj|2008h . The model 
by De Lucia & BlaizoT i 2007 ) predicts a constant merger rate evo- 
lution, visibly at odds with the observations and the Bower et al. 
model results. 

In Fig. [3] we investigate the dependence of the merger rate 
on the mass ratio of the pair components. We show results for the 
Bower et al. model from which we selected galaxy pairs with mass 
ratios smaller than 4:1, 6:1, and 10:1, as described in Section [3] 
The amplitude of the merger rate derived from the model depends 
strongly on the mass ratio adopted to select the pairs, and it is com- 
patible with the observed one only when we consider a mass ratio 
threshold of 10:1. The fit of the results with the same function as 
in equation l|2]l gives m — 0.59 for merger rates derived from pairs 
with mass ratios smaller than 10:1, 0.63 for mass ratios of 6:1, 
and 0.68 for mass ratios of 4:1. These results are consistent with a 
roughly flat evolution of the merger rate independently of the mass 
ratio threshold adopted. 



Figure 3. Evolution of the galaxy merger rate predic ted by the semi- 
analytical galaxy formation model of Bower et al. (2006) considering mass 
ratios of 4:1, 6:1, and 10:1 to select galaxy pairs, as described in Section[3] 



5 SUMMARY AND CONCLUSIONS 

In this Letter, we have explored galaxy mergers in model universes 
drawn from semi-analytical galaxy formation models implemented 
on the Millennium Simulation of structure growth. Fr om the public 
alaxy cat alogues and merger t rees g enerated by the IBower et al.1 
2006) an d l5e Lucia & Biaizod Emm models we retrieved the pair 
fractions and the galaxy merger rate (number of galaxy mergers 
per Gyr), and compared their redshift evolution with observational 
results in the redshift range < 2 < 2. 

The results from our study indicate that the evolution of the 
galaxy merger rate has been discrete since z ~ 2, following an 



exponential function given by T„ 



(1 



with m rang- 



ing from 0.6 to 0.8 depending on the luminosity and mass ratios 
of the merging galaxies. We find that this be haviour is true only 
for the Bower et al. model, since the model bv lDe Lucia & Blaizotl 
d2007h predicts a null merger rate evolution (m ~ 0), at odds 
wit h the observations and the Bower et al. trends. The results for 
the |Pe Lucia & Blaizot 1 20071) mo del can be related to the study 
carried out bv lGuo& White! d2008h on the galaxy growth through 
mergers and star formation based on it, where they have found a 
nearly constant growth by mergers in the range < z < 2, inde- 
pendent of the galaxy stellar mass. 

The weak evolution of the galaxy merger rate found in our 
analysis is consistent with recent observational results obtained by 
different approaches which give supp ort for a flat evolution of the 
galaxy merger rate, wi th m ~ 0.5—1 l lBundv et al.l2004l ; lLin et all 
l2004l : lLotz et al.l2008l) . Previous studies have shown that in the hi- 
erarchical framework the merger rate of cold dark matter haloes 
increases rapidly with reds hift also as (1 + z) m , but with the 
exponent m ~ 2 — 3 (e.g. iGott lober. Kly pin & Kravtsov]|200ll: 
iFakhouri & Mai |2007|) . However, as discussed by iBerrier et al .1 
d2006|) . multiple galaxies may occupy the same parent dark matter 
halo, resulting that the predicted galaxy merger rate should evolve 
slower than the halo merger rate. Actually, in the present work we 
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have confirmed this tendency. The weak evolution of the galaxy 
merger rate predicted by the galaxy formation models shows that 
the mass assembly evolution of galaxies through mergers does not 
follow the rapid evolution of the halo merger rate obtained in pre- 
vious studies. 
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